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Abstract 

We analyze the stability properties of equilibrium solutions and periodicity of orbits in a two- 
dimensional dynamical system whose orbits mimic the evolution of the price of an asset and the 
excess demand for that asset. The construction of the system is grounded upon a heterogeneous 
interacting agent model for a single risky asset market. An advantage of this construction 
procedure is that the resulting dynamical system becomes a macroscopic market model which 
mirrors the market quantities and qualities that would typically be taken into account solely at 
the microscopic level of modeling. The system's parameters correspond to: (a) the proportion 
of speculators in a market; (b) the traders' speculative trend; (c) the degree of heterogeneity 
of idiosyncratic evaluations of the market agents with respect to the asset's fundamental value; 
and (d) the strength of the feedback of the population excess demand on the asset price update 
increment. This correspondence allows us to employ our results in order to infer plausible causes 
for the emergence of price and demand fluctuations in a real asset market. 

The employment of dynamical systems for studying evolution of stochastic models of socio- 
economic phenomena is quite usual in the area of heterogeneous interacting agent models. How- 
ever, in the vast majority of the cases present in the literature, these dynamical systems are one- 
dimensional. Our work is among the few in the area that construct and study two-dimensional 
dynamical systems and apply them for explanation of socio-economic phenomena. 
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1 Introduction 



We investigate equilibria and stability properties of the discrete time dynamical system 
generated by planar map in the following manner 

/ \ / \ / Pn-l + A(i„_i 

^"^"-^ ^"^^-'^ \ +(l-a)[l-2$(p„_i + A4_0] 

In the dynamical system ([1]), a G (0, 1), A and J are positive real numbers and $ is a 
probability distribution function. Since we do not impose rigid constraints on this dis- 
tribution function, the analysis of ([T]) becomes a nontrivial task. Nevertheless, our study 
of ([1]) was motivated by its potential applications rather than by intrinsic mathematical 
challenges yielded by the non-linearity. The applications stem from the link of the sys- 
tem to a stylized model of a single risky asset market represented by a stochastic process 
constructed by us. 

The stochastic process just mentioned is a microscopic model of a market in the sense 
that it mirrors individual behavior of each market agent. The model will be explained 
in Section 13. ![ We constructed it using ideas from the area of Heterogeneous Interacting 
Agent Modelqj (HIAMs). HIAMs have a long history, wide applications and employ meth- 
ods from a variety of disciplines including Statistical Mechanics and Interacting Particle 
Systems (see f?]). The experience showqj that HIAMs are able to adequately describe 
macroscopic characteristics of socio-economic systems by means of mirroring system's mi- 
croscopic components and their mutual interactions. Our HIAM was constructed with the 
aim to investigate evolution of two particular macroscopic characteristics: the asset price 
and the population excess demand for an asset in real world markets. In the constructed 
model, the evolutions of interest are represented by stochastic processes. These processes 
possess the following property: as the model's agent population increases, their trajecto- 
ries converge to orbits of the dynamical system ([1]) (this convergence will be explained in 
Section [3^ . This convergence yields the following interpretation of the parameters and 
variables of ([1]) in terms of a single risky asset market: 

Pn corresponds to the price of the market asset at time n; 
dn corresponds to the population excess demand for the asset at time n 
(in the sequel, it will be called excess demand for short); 

a corresponds to the proportion of speculators among the market traders; 

J corresponds to the traders' speculative trend; 

A corresponds to the feedback of the excess demand on the increment 
of the update of the asset price; 



^The works [T] and [4] provide an broad survey of the history of the HIAMs, as well as of the current 
statc-of-thc-art in the area. In particular, our HIAM shares many features in common with the HIAM 
from the seminal works of Kirman ([5]) and of Lux ([8]). Among the more recent studies that present 
HIAMs similar to ours, we wish to single out [2] and [9] because the problems addressed by them are 
related to the issues discussed by us here. However, we shall not pursue a detailed comparison of our 
model to the already existing ones since this is not essential for our presentation. 

^See the review works U and [1], and citations therein. 



2 



$ corresponds to the distribution of the deviations around a constant v 

of individual evaluations of the asset fundamental value, where v corresponds 
to the market fundamental value of the asset. 

This interpretation will be explained in details and justified in Sections 13.11 and 13.21 It 
will allow us to re-phrase our results so that they explain whether and why an asset price 
and excess demand in an asset market would or would not oscillate as time goes on; this 
will be the contents of Section 13.31 

The above mentioned convergence of HIAM's trajectories to orbits of an appropriate 
dynamical system is a well known fact in the HIAM area. It is usually employed, § but 
almost always in such a way that the resulting dynamical system is of the dimension 
one. This is motivated by the wish to use the graphical analysis as a tool for study 
this system. Contrasting, in the present work, the limiting dynamical system is two- 
dimensional. The dimension increase happens to make a difference because it allows one 
to see socio-economic phenomena that are invisible through the lenses of one-dimensional 
dynamical systems. The price to pay for this is the growing complexity of the proofs. 

The results of our study of ([T]) are formulated in Theorems [T], [3] and H] of Section [2l 
Theorem [1] asserts that the point (0, 0) is the unique equilibrium state of the dynamical 
system ([T]) and determines when it is locally asymptotically stable and when it is unstable. 
An interesting feature of this theorem is that it characterizes the stability/instability via 
a relation between just two expressions involving the model's parameters (that are, we 
recall, three real numbers and a probability distribution function). This reduction of 
the parameters' space allows us to draw a phase diagram (Figure [T]) representing the 
stability/instability of the unique equilibrium point. Our second result. Theorem [21 gives 
conditions on the parameter values for the appearance of a stable periodic orbit of the 
dynamical system. Generalizations of this theorem are indicated in Conjectures [1] and [21 
Our third result. Theorem [H identifies parameter values for which the point (0, 0) is a 
globally stable equilibrium. 



2 Results 

Here, we present results of our study of the dynamical system ([T]) in which 

a E (0,1), J > 0, A > 0, (2) 
and $ is a probability distribution function that satisfies the following conditions: 



(a) $ is everywhere differentiable; we shall denote its derivative by 

(b) sup,{$'(a;)} = $'(0); 

(c) $' is strictly increasing on (— cxo, 0), and strictly decreasing on (0, oo); 

(d) $(0) = l/2; 

(e) a random variable with the distribution $ has zero expectation; 
(/) a random variable with the distribution $ has finite variance. 



(3) 



■^See [6] for a general account on the employment of this convergence for revealing socio-economic 
phenomena. 
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We might have relaxed significantly the constraints ([3]) on $, but this would not have 
broadened the application of our results and would make their proofs more cumbersome 
without introducing essentially new ideas. We note also that for applications, it is rea- 
sonable to admit that $ is a Gaussian zero mean distribution function. This function 
satisfies the above conditions. 

2.1 The uniqueness of equilibrium and its basic stability prop- 
erties 

Theorem 1 (uniqueness of equilibrium and its stability). 

(a) The origin (i.e., the point (0,0)) is the unique equilibrium of the dynamical system 

(b) If X^'{0) > 1 + 2q;J$'(0)), then the origin is an unstable equilibrium of (Qp, while 

A$'(0) < 1 + 2q;J$'(0), then the stability of the origin depends on the value of 
2a J$'(0) in the following manner: 

if 2aJ(^' (fS) < 1 then the origin is locally asymptotically stable; 
z/2q;J$'(0) > 1 then the origin is unstable. 

Proof of ThmUl If {pn,dn) = {po,do), for n E N, then in virtue of the first equation of 
([1]), Xdo = 0. This implies that do = because A > by construction. On substituting 
do = in the second equation of ([T]), we get that 1 — 2$(po) = 0. This can hold only if 
Po = because $ is monotone and $(0) = 1/2 (the properties ensured by the assumptions 
([3])). This completes the proof of item (a). 

We proceed with the proof of (b) of the theorem. From ([1]) we easily get the Jacobian 
matrix of at the origin: 

( ' ' 

\ -2$'(0) 2(aJ- A)$'(0) 

Its characteristic polynomial is /i^ — [l + 2(aJ — A)$'(0)]/i + 2aJ$'(0) whose roots are 

/ii = i[l + 2(aJ-A)<l>'(0) + yA], 
/i2 = |[l + 2(aJ-A)<l>'(0)-yA], 

where 

A := [1 + 2{aJ - A)$'(0)]^ - 8aJ<l>'(0). (6) 

The assertions in (b) will all follow from the principle of linearized stability via analysis 
of the values of fii and fi2- The analysis are split into six cases. They may be seen from 
Fig. [TJ This figure is a kind two dimensional phase diagram because - due to the form 
of the roots - the stability analysis relies on the relation between just two expressions: 
2aJ<l>'(0) and 2A$'(0). 
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w = 2A$'(0) 



Region 6 




[^initial) ^initial y 

.Region 3 
Region la 



(Wfinal, Wfinal) 



Region lb 



V Region 2 




Region 4 



u = 2aJ$'(0) 



Figure 1: Phase diagram representing the stabihty of the unique equihbrium of the dy- 
namical system ([1]). 



AU-throughout below, we shall usually use in our calculations that A > (the con- 
straint on A imposed by our construction) and that $'(0) > (the inequality that follows 
from ([3]-a, b, c)). 



Case 1: A < (region 1 in Fig. [T]). 

We have: 



l/iil' = |/i2|' = ^(1 + 2{aJ - A)<l>'(0))2 + (v^)' = 2«J$'(0), 



(7) 



and therefore the origin is locally asymptotically stable, in case 2q;J$'(0) < 1 (region la) 
and is unstable, in case 2aJ$'(0) > 1 (region lb). 

Case 2: A > 0, 2A$'(0) < 1 + 2aJ$'(0) and 2aJ$'(0) < 1 (region 2 in Fig. [I]). 
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From the hypotheses, we obtain immediately: < l + 2(aJ — A)$'(0)) < 1 + 2q;J$'(0). 
By adding and subtracting (l + 2a;J$'(0)) , we get: 

f 

Thus, < 1 - 2aJ<l>'(0) and then: 



A = (l + 2aJ<l>'(0))^-8aJ<l>'(0)+ [l + 2(aJ-A)<l>'(0)]^- (l + 2aJ$'(0))' 



1 + 2(aJ- A)$'(0) (l-2aJ$'(0)) 

< 2 + 2 

= 1-A<I>'(0) 

< 1. 

Since the coefficient of fi in the characteristic polynomial is negative we conclude that 
1/^2 1 < lAiil < 1 and hence the origin is locally asymptotically stable. 

Case 3: A > 0, 1 + 2a J$'(0) < 2A<I>'(0) < 2 + 4a J<l>'(0) and 2a J<l>'(0) < 1 (region 3 in 
Fig. [ID . 

By simple algebraic manipulations, it follows that —(1 + 2aJ$'(0)) < (1 + 2(aJ — 
A)$'(0)) < 0. Therefore, as above, we obtain: a/A < 1 — 2aJ$'(0) and consequently 

1 + 2(aJ- A)<l>'(0) (l-2aJ<l>'(0)) 



1^-2 > 



2 

4aJ$'(0) - 2A<I>'(0) 



> -1. 

Since the coefficient of n in the characteristic polynomial now is positive we conclude that 
\fJ'i\ < \f^2\ < 1 and hence the origin is locally asymptotically stable. 

Case 4: A > 0, 2A$'(0) < 2aM'{0) - 1 and 2aJ$'(0) > 1 (region 4 in Fig. ^. 
Now, we have 2(a J — A)$'(0) > 1 and, therefore: 

l + 2(aJ-A)1>X0) , r- 
Aii > hvA 

1 2(aJ-A)1>X0) 

2 2 

> 1, 

proving that the origin is unstable in this case. 

Case 5: A > 0, 1 + 2a J$'(0) < 2A$'(0) < 2 + 4aJ$'(0) and 2aJ<l''(0) > 1 (region 5 in 
Fig. [ID . 

Since A > 0, we must have either 2A$'(0) > (l + ^/2aJ^'{0)Y or 2A<I>'(0) < (l - 

^2aJ$'(0))l But, since (l - y^2^M^y < 1 + 2aM'{0) < 2A$'(0) then the second 
alternative is false. Thus: 

2A$'(0) > (l + v/2aJ$^)^ 

= l + 2^/2^M^ + 2aJ^\0) 
> 3 + 2aJ$'(0). 
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We then obtain 2[aJ — A)$'(0) < —3, from which it follows that 

1 + 2(aJ- A)<l>'(0) 

< ^ 

< -1, 

proving instability. 

Case 6: A > 0, 2A$'(0) > 2 + 4aJ<l>'(0) (region 6 in Fig. [I]). 

In this case, since l + 2(a J- A)<l>'(0) < l + 2aJ$'(0)-2-4«J$'(0) = -l-2aJ$'(0) < 
0, we have A = (1 + 2{aJ - \)^'{())f - 8a J > (1 + 2aJ^'{Q)f - 8aJ = (1 - 2a J$'(0))2. 

If 2aJ^'{0) < 1, then > 1 - 2aJ<l>'(0), and we obtain 

l + 2{aJ - X)^'{0) 1 - 2aJ^'{0) 

2 2 
2aJ$'(0) - A$'(0) 

2aJ$'(0) - 1 - 2aM'{0) 

-1. 

If 2aJ$'(0) > 1 then VA > 2aJ<l>'(0) - 1, and we obtain 

1 + 2(aJ- A)'l'XO) _ -l+2aM'{0) 
^ 2 2 

= 1-A<I>'(0) 

< 1 - 1 - 2aJ<l>'(0) 

< -1. 

Thus, /i2 < 1 in both case, and the instability follows. □ 
2.2 Hopf bifurcation and periodic orbits 

The occurrence of Hopf bifurcation in the dynamical system ([T]) is the issue of Theorem [2] 
and Conjectures[T]and[2]of the present section. All they are based on Theorem [2] below. It 
contains a description of the Hopf bifurcation phenomenon and provides conditions that 
are sufficient for it to occur. This bifurcation implies the emergence of periodic orbits in 
the dynamical system for particular set of its parameters' values. 

Theorem 2 ([3], page 474, Poincare-Hopf-Andronov theorem for maps). 
Let 

F : Rx (Rf ^ R^; (r/, x) ^ F(r/, x) 
be a map depending on a real parameter rj satisfying the following conditions: 

(i) F{ri, 0) = for rj near some fixed rjo; 

(a) DF('r],0) (that is, D,^F{r],x)\^_^) has two non-real eigenvalues n{r]) and fi{r]) for r] 
near rjQ, with |/i(77o)| = 1; 



< 
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(iv) /i'=(r/o) 7^1 /or A; = 1,2, 3, 4. 

Then there is a smooth rj-dependent change of coordinates bringing F into the form 

F(r7,x)=^(r7,x)+0(||x||5) 

and there are smooth functions A{ri), Birf) and ^l{r]) so that in polar coordinates the func- 
tion J-'{ri, x) is given by 

( r\ ( \ii{vi)\r - A{r])r^ \ 

\e J \e + Q{7]) + B{7]y J ■ 

If A{riQ) > 0, then there is a neighborhood U of the origin and a S > such that, for 
\v — Vol < S o-nd X G t/, the uj— limit set of ^ is the origin if f] < rjo, and belongs to a 
closed invariant curve T{r]) encircling the origin if rj > rjQ. Furthermore Tirio) = 00 

If A{r]o) < 0, then there is a neighborhood U of the origin and a 6 > such that, for 
\v — Vol < ^ o-nd X G f/, the uj— limit set of ^ is the origin if v > Vo, <md belongs to a 
closed invariant curve T{ri) encircling the origin if v < Vo- Furthermore r(?7o) = 0. 

Remark 1. We present here the method that we shall employ to calculate ^(770)- The 
presentation follows [3] . If the linear part of map F at ?7o is written in the Jordan canonical 
form, 

where a and h relate to the eigenvalue ^{rj) via the equation ^{vo) = a + ib, then for A{r]) 
defined in ([H]), it holds that 

+ ^ kill + 1^02!- Re (10) 

where (all the derivatives below are evaluated at (xi,X2) = (0,0)) 

C20 = I {{Gl)xixi — {Gl)x2X2 + 2(G'2)xiX2 + i[{G2)xiXi — {G2)x2X2 — '^iGi)xj^x2]} ! 
Cll = i {{Gl)xixi + {Gl)x2X2 + i[{G2)xixi + (^2)^2x2]} ; 

^02 = I {{Gl)xixi — {Gl)x2X2 — '2{G2)xiX2 + ^[(G'2)xixi " {G2)x2X2 ~ '^{Gl)xiX2]} ! 

and 

C2I = Yq {(^1)2:12:12:1 + {Gl)xiX2X2 + {G2)xiXiX2 + {G2)x2X2X2 

+ i[{G2)xiXiXi + {G2)xiX2X2 — {Gi)xixiX2 — (^1)2:22:2x2]} • (11) 

"'The phenomenon described in the last two sentences is called Hopf bifurcation. We add to this name 
the term supercritical in order to distinguish it from the case described in the next two sentences that 
will be called the suhcritical case. 



A{7]o) = Re 



1 - Kvo) 



20 
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Theorem [3] below reveals the occurrence of Hopf bifurcation in the dynamical system 
([T]) when it satisfies certain additional assumptions. How general these assumptions are 
will be discussed after the theorem's proof. The same discussion will expose at an intuitive 
lever the ideas behind the proof. 

Theorem 3. 

Suppose the dynamical system ^ satisfies the following assumptions: 

Assumption 1. The parameters X,a,J and the distribution function $(•) are all 
functions of a single real variable rj that are defined on some (nonempty) inter- 
val [?7initiai, VRnai]; thc functions Ar,,ar), Jr] and the distribution functions are 
smooth enough to ensure that \E', the map that generates the system is in 
{pn-i,dn-i) and t]. 

Assumption 2. The functions a-r^, J. an(i$^(0) (the latter means the value of d^r]{x)/dx 
at X = 0) are all increasing in t] on the domain [?7initiai, ''^finai]- 

Assumption 3. (l + 2(a^J^-A^)<l>'^(0))2-8a,,J,,<l>'^(0)) < 0, for each rj E [r^mitiai, ??finai]- 

Assumption 4. There exists a unique rjo E [T^mitiai, ''7finai] such that 2q;^(, J^(,$^^(0) = 1, 
and r]o ^ ij^^itiai, Vo 7^ '7finai- 

Suppose in addition that 

ia)K = l^ (6)$;(0) = 1, (c)$';(0) = 0, (rf)$:;;,(0)<0. (12) 
Then the system undergoes the supercritical Hopf bifurcation when rj passes through 

Vo- 

Proof. In virtue of Thm. [2] and Remark [H in order to prove the present theorem, it is 
sufficient to show that the mapping that generates the dynamical system ([1]) satisfies 
the conditions (i) - (iv) of Thm. [2] and that v4(?7o) is a positive number. We do so in 
Steps 1-5 below. 

Step 1: Since Theorem [1] ensures that (0, 0) is an equilibrium point of ([T]) for any values 
of the parameters of this dynamical system, then \l/ satisfies (i) for each 77 G [77initiah ''^finai]- 

Step 2: We start recalling facts and results from the proof of Thm. [1] that we shall 
need below: (a) the Jacobian matrix of the mapping \E' at (0,0) (i.e., DF{ri,0), in the 
notations of Thm. [2]) was calculated and the result is presented in (b) the matrix' 
eigenvalues were calculated and their expressions are presented in ([2]); (c) it was proved 
that these eigenvalues are non-real numbers provided A < 0, with A being defined by 
(Q; (d) the modulus of each eigenvalue was found to be equal to ^/2xyJW(0), provided 
A < 0. 

Now, since Ass. 3 ensures that A < for every rj, then (c) above implies that 
satisfies the first part of the condition (ii). As for the second part of this condition, it is 
implied by Assumptions Sand 4 and the the fact (d). 

Step 3: The fact (d) from the list of Step 2 and Assumptions 2 and 3 imply straight- 
forwardly that the condition (iii) is satisfied by the eigenvalues of \E' for every rj, and in 
particular, for r/o- 
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Step 4- From the fact (b) of the hst of Step 2 and from Ass. 3, it follows that 

/^(^o) = ^ (13) 

J8a,„^$;„(0)) - (1 + 2K„^ - 
+^ 2 • 

Applying Ass. 4 and constraints (fT2|) to (fT3i) . we get that /i(77o) = i + From this, 
H^ir^o) ^ 1, k = 1,2, 3, 4, proving hence the validity of condition (iv). 

Step 5: In this step, we shall prove that A{riQ) > 0. In the calculations that follow, all 
functions that depend on t] will be evaluated at tjq. This allows us to omit r]Q in the nota- 
tions throughout the proof, namely, we shall write A, X,a, J and $'(0) for, respectively, 
A{r]o), Xrjo,arjo,Jrjo and $^(,(0). 

Since we intend to use the expression ffTOl) then we need to get a Q-like form of the 
map We shall get it from the Jacobian matrix at (0, 0) of the map that we've 
calculated in the proof of Thm. [T] (the matrix is presented in (jlj)). Since Ass. 4 imposes 
that 2aJ^'{0) = 1 then this matrix acquires (recall, at r^o) the following form: 



M 



1 A 

-2<l>'(0) 1-2A$'(0) 



It is then easy to check that the matrix 



1 
-$'(0) 



P={ ^un. v/T^ )' With £:=1-A$'(0), 



puts M in the Jordan canonical form: P ^MP = ( ^ ^ ) . Therefore 



in the new variables ( ^ ^ ^" , n G N, the map \I' acquires a Q-like form: 

f ] = P-'MP ( ) + P-' (^ - M) P . , 



e 



- \ I Pn-1 \ , I Gi{pn~l, dn-l) 



where 

Gi(xi,X2) = 0, 

— £ 

+$'(0) (3 - 2a<l>'(0)) xi - X-^ (1 - 2a$'(0)) Vl -e^xa}, 

and where 



i^i = 1 - (A - J)$'(0), K2 = A-Vrpi2(A - J), 

K3 = 1-A$'(0), K4 = VT^. 
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From the expressions for Gi and G2, we shall now obtain i^^^'s following the formulas 
from Remark [H Differentiating G2 we get : 



a2 2\ 

—G2{xi,X2) = ; {aKf^"{KiXi + K2X2) + (1 - a)Kl^"{K^xi + KiX2)) . 



dx\ ^' v^r^ 

From this result, in virtue of the assumption flT^c). we conclude that ^6*2(0,0) = 0. 
By similar arguments, we conclude that all second order derivatives of G2 vanish at 
(a;i, X2) = (0, 0). These facts and the fact that Gi = imply that ,^20 = ^11 = ^02 = 0. In 
order to find ^21, we calculate third order derivatives of G2- 

6-2(0,0) = - [aKf + (1 - a)Kl) 



dXn 



2W" 


(0) 


VI- 






(0) 


VI- 




2A<I>"' 


(0) 


VI- 




2A<I>'" 


(0) 


Vi- 


£2 



dxf 

^^G2(0,0) = -^=maK',K2 + il-a)K!K,) 

^^G2(0,0) = -^=maK,K! + il-a)K,K!) 

G2(0,0) = =^{aK| + {l-a)Kl) 



Now, from the relations 2aJ$'(0) = 1 (that holds because of Ass. 4) and $'(0) = 1 
(ensured by (fT2l-b)) we get that a = (2J)^^ This allows us to substitute a by (2J)~^ in 
the expressions for the derivatives of G'2. We also substitute there A and $'(0) by | and 
1, respectively (these substitutions are justified by (IT^a. b)). The resulting simplified 
expressions for the derivatives lead, via the formula (ITT]) , to the following: 

^21 = |^[(-20J2 + (13v^-4)J + 3v^-10) 
+i(-10 J2 + (4v^ _ 7) J + 2^3 - 5)]. 

Plugging in flTU]) the expressions for .^^^'s obtained above, we finally get that 
A = -Re(/i6i) = -$'"(0)^±-^[5j2 + (1 - 2V3)J + 1]. 

This expression for A and the assumption (fT2] -d) ensure that A > for every J. This 
completes Step 5 and the proof of Theorem [31 □ 

We proceed with the discussion of how general the Hopf bifurcation phenomenon is 
for the dynamical system ([1]). 

In order to discuss Hopf bifurcation in a dynamical system with the help of Theorem [2l 
the minimal necessary condition is that all system's parameters be expressed as functions 
of a unique variable. In the framework of Theorem[31 this condition is ensured by Assump- 
tion 1. Let us accept it now and let us consider then two functions: := 2a^J^$'^(0) 
and := 2A^$'^(0). We recall from the proof of Thm. [T]that (a) the quantity denoted 
there by A can be expresses as a function of solely u and w, (b) when A{u,w) < then 
the eigenvalues of the linearization of \I' at its fixed point (0, 0) are non-real numbers. 
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and (c) when u = 1 then the modulus of each eigenvalue is 1. The facts (a,b,c,) ensure 
that if the parameter t] introduced in Assumption 1 is such that A(m^, ly^) < 0, for all 
7], and if = 1, for some rjo, then the property (ii) of Thm. [2] is satisfied. These "if" 
conditions are provided by Assumptions 3 and 4. As for the Ass. 2 ensures the validity 
of (iii) of Thm. O For this, Assumption 2 may be not the minimal sufficient condition, 
but we did not search for such. Turning our attention to the condition (i) of Thm. [21 we 
easily see that it is valid for \E' due to ThmIT](a) without the necessity for any additional 
assumptions. 

We thus have showed that Assumptions 1- 4 ensure the validity of conditions (i) - 
(iii) of Thm. [2l It seems to us that these assumptions are also sufficient for the validity 
of the condition (iv) of Thm. [2] and the inequality A{rjQ) ^ 0. Accordingly, we formulate 
the following: 

Conjecture 1. 

The dynamical system (J\) exhibits periodic orbits, when the initial point is sufficiently 
close to (0, 0) and when the parameters are such that the corresponding point {u, w) is 
sufficiently close to the right of the interval {(u,w) : u = 1, I < w < 4}|§ precisely to 
state, there is a neighborhood U of (0, 0) such that the uj-limit of any orbit of (QP with 
initial condition in U belongs to a closed curve encircling (0,0), provided 2aJ$'(0) is 
slightly bigger than 1 and A$'(0) is between and 2. 

We did not pursue proofs in the degree of generality that would allow us to justify 
rigorously the generic property formulated in this conjecture. Rather, we considered 
two particular cases. The first case is presented in Thm. [31 There, we assumed an 
additional constraint f[T^ that helped a lot to simplify the calculations needed to establish 
the condition (iv) and the inequality A{rio) > 0. The second case is presented in the 
statement below. This case attracted our attention because it arises in applications of 
our mathematical study of the dynamical system The reason for this is explained in 
Section [331 

Conjecture 2. 

// the parameters of the dynamical system (QP satisfy Assumptions 1-4 '^^^ ^ Normal 
Distribution with null mean then the system undergoes the supercritical Hopf bifurcation 
when rj passes increasingly through tjq. 

In the case of Conjecture [21 the analytic verification of the condition (iv) and the 
inequahty A{rio) ^ turned to be an extremely tedious task. To carry out this task, we 
resorted to numeric methods: we verified the condition and the inequality numerically for 
a grid of parameter values. 

'"'If the point is close to this interval on its left, then (0,0) is locally asymptotically stable equilibrium 
of the dynamical system (|T]) and hence the periodic orbits described in this statement cannot occur. 
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2.3 Global asymptotic stability for small values of the parameter 

A 



Theorem 4 (a sufficient condition for global asymptotic stability). 
Suppose 2aJ$'(0) < 1. Then, there exists a positive number Ac such that if \ E (0, AJ 
then the equilibrium (0, 0) of the dynamical system (QP is globally asymptotically stable. 
(The numeric value of Ac will be specified in the proof by [JW-) 

Remark 2. On comparing the assertions of Thms. [1] and H] one concludes that the 
assumptions of Thm. |4] must be a sub-case of the assumptions of Thm. [H This inclusion 
can be established by simple calculations grounded on (HOI-fii)). one of the conditions that 
will determine the value of Ac- We shall omit these calculations. 

We now proceed with the argument that will lead to the proof of Thm. HI The 
arguments employ essentially the assumption ([2]) that ensures that J > and a > and 
thus, allows us to divide by J and a. Actually, the theorem's assertion remains true even 
when either J and a is equal to 0, but this case requires a specific argument which will 
not be presented here. 

Upon the substitution of Pn-i + A(i„_i by Pn in the second equation of the dynamical 
system ([1]), it acquires the form that makes it clear that the passage from 
to {pn,dn) can be considered as a two-step procedure: first, pn is obtained from Pn^i 
and dn-i via the first equation of the system, and then, dn is obtained from dn-i and 
Pn via the second equation. The second equation, or better to say, its non-linearity is 
the principal cause for the difficulty of studying the evolution of the dynamical system 
([1]). In the present proof, we overcome the difficulty by getting convenient estimates of 
two quantities (to be defined in fl23|) and flM|) below) that relate to this equation and 
characterize its features that are important for our analysis. The presentation of these 
estimates becomes more transparent - in our opinion -, when we express the equation 
through an appropriate one-parameter M ^ M mapping in which p„ is the parameter 
and dn-i is the argument. In order to do so, for the "parameter" p G M, we define 
9pi') '■ (—00,00) — i> [—1, 1] as follows]^ 

^p(rf) =a[l-2<l>(p- Jrf)] + (l-a)[l-2$(p)], rfG (-00,00). (14) 

Now, we can re-write ([1]) in the desired form: 

Pn = Pn-l + A(i„_i (15) 
dn = 9p„ (dn-i) (16) 

From now on, we shall use the shorthand notation: 

b := 2aJ$'(0); (17) 

note that 6 > (because of ([3]-(a), (c))) and b < 1 (because of the theorem's assumption). 

^Although the argument of gp{-) cannot exceed 1 in modulus - since it has been defined as the 
mathematical analogue of population excess demand -, it turns out to be convenient to extend the 
function's argument domain to the whole M. 
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We start with a list of basic properties of gp{-). For each p, 



o (y(p(-) is everywhere differentiable function and g'pid) = 2aJ$'(p — Jd); (18) 
o sup{^?;(rf)} = (4) = 2aJ<l>'(0) = 6 G (0, 1); (19) 

o in particular, as a consequence of (ITSl) and (IT^ . 5'p(-) is a contraction: 

kp(a;) -^p(l/)| < < (20) 

o Qpi') a monotone strictly increasing function; (21) 
o the graph oi gp{-) passes through the point 

(rf*, y*) with d* = ^ and y* = {I - a) [1 - 2<l>(p)] (22) 

that satisfy / ^ ^ (>) ^ ^ > (>) 0, y* < (<) 0, 

that satisfy. | ^ < (<) q ^ \ < (<) 0, y* > (>) 0. 

All these properties stem from f|T^ in combination with the properties of $ assumed in 
([3]); exactly to state, we used (^a) for ([IHD, ([lb) for ([19]), ([Ic) for ([2l]), and ([Id) for 



We now define the quantities V^p) and 7^(p) that will play central role in our argu- 
ments: 

V^p) := the solution of gp{d) = d, hence, gp(V{p)) = Vij)), (23) 

{— oo, when gp{d) > 0, Wd, 
+00, when gp{d) < 0, Wd, (24) 
the solution of gp{d) = 0, otherwise, and hence, gp{JZ{p)) = 0. 

Figure [2] helps the visualization of these definitions. We observe that the properties f[T9l) 
and fl2T|) ensure that T'(p) and 7^(p) are well and uniquely defined. 
We shall frequently use the following properties of 'D{p) and TZ{p): 

o if p > then 7^(p) > - > 0, and if p < then 7^(p) < - < 0; (25) 

J J 

o iip<p' then V{p) > V{p'); (26) 
o if p > then V{p) < 0, and if p < then V{p) > 0. (27) 

(!25|) follows directly from (I2T!) and (!22|) . To prove (!26l) . we argue as follows: First, from 
the monotonicity of $ and the definition ([HI), we get: 

o gp(-) is decreasing in p, i.e., gp{d) > gp/{d)Wd, when p < p' . (28) 

Then, from this property and the fact that |5'p((i)| < 1, we conclude that the first coor- 
dinate of the intersection point of the 45° line in the plane with gp{-) decreases with the 
increase of p. This conclusion is exactly the property fl2B]) . As for (1271) . it follows from 
fl2^ because 1^(0) = (this equality holds because (7o(") passes through (0,0) as ensured 
by (^d) and ([IID). 

There are two more properties of T>{p) and 7^(p) that we shall need for the proof of 
Theorem |H These are provided by Lemmas [1] and [2] below. 
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Figure 2: Illustration to the definition of the quantities T>{p) and TZ{p) and the arguments 
that employ these quantities and their properties. 

Lemma 1. If p > then V{p) > ~j^^p, and if p < then V{p) < -j^^^p. 

Proof. We prove the case p > 0; in the case p < 0, the proof is similar. 
We shall need the following estimate: 

1 - 2*(p) l^-S'™ -2 f\'(.)d. -2*'(0)p'S' (29) 

Jo OiJ 

No-ftQ, we construct the tangent line to gp{-) at the point {d*, y*). From (IT5I) . (IT^ and 
(122]) . we have that the slope of the line is b and the graph of gp{-) to the left of {d*, y*) lies 
above this line. We also have that d* > and y* < (because of fl22l) and the assumption 
p > 0). Hence, 'Dip), the coordinate of the intersection of gp{-) with the diagonal, is 
larger than T'(p), the coordinate of the intersection of the tangent line with the diagonal. 
By simple calculations, we get that V^p) = (1 — b)~^ |(1 — a) [l — 2$(p)] — 
Combining this with the inequality 'P(p) > T^lp) just obtained and with the estimate 
fl2U]) . we complete the proof: V{p) >V{p)> ^ = -j^^^p. □ 



''The following constructions are illustrated in Figure [2l 
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Lemma 2. Let (a defined below might have been any other number from (0, b)): 



a := 6/2. (30) 

Let next Xe < 0, Xr > be such that: 

2aJ^'{x) > a, Wx E [xe,Xr] (31) 

(the existence of xi and Xr follows from the properties (O-a, b, c) of ^; in fact, xi and Xr 
are the negative and the positive solutions of the equation 2aJ$'(x) = a). Finally: 

h := mm — , lj> x min Xr-}. (32) 

Then, TZ{-) satisfies the following properties: 
(a) Ifp > then TZ{p) > ^fp, and if p < then n{p) < -^fp. 
{h) If{]<p<h then n{p) < ^^p, and if-h<p<0 thenU^p) > -^^p. 

Proof. We shall prove both (a) and (b) for p > 0. For p < 0, the proofs are analogous. 

Due to (132|) . a < p < h then x G [— 0] implies x G [a;^,0]. This implication and 
( I3T1) ensure that $'(x) > (2aJ)~^a, in case x G [—p, 0] and p G [0, h]. We use this fact to 
derive the following inequality: 

1 - 2$(p) ^" ^=^' ^''^ -2 r ^'{x)dx< - —p, for all p G [0, h]. (33) 

Jo 

Recall now from fl22]) that the graph of gp passes trough the point {d*, y*). We construct 
two straight lines passing through this point and having slopes b and a, and we denote 
by T^iower and 7?.upper the coordinates of their respective intersections with the d axis 
(see Figure [2]). Directly from this construction, we have that 7?.iower = d* — ^ and 
7^uppcr = d*-^. Since y* = (l - a)[l - 2^p)] then: 



p ]^ — cx Qj (jooji p ]^ — (y h 

Slower > ^ H ^tP, Vp > 0, and 7^uppcr < ^ H 7 — P, Vp G [0, h]. (34) 

J aJ b J aJ a 

Next, due to ([32D, if < p < then x G [0, ^^p] implies that a; G [0, x,]. This 
imphcation and (l3TI) ensure that 2a J$'(x) > a, in case x G [0, ^^^p] and p G [0, h]. This 
conclusion and the relation ( fTSi) between (yf^ and $' together ensure that if p G [0, h] then 
fl'p('^) > a for each ci G [j, j + ^^^p]- From this inequality, it follows that the graph 
of gp{d), G [j, J + ^^^p]; liss above the hne of slope a constructed by us, provided 
p G [0, h]. This fact implies that 7l{p) < T^upper, when p G [0, h] (for this implication to 
be valid it is important that the interval on which gp{-) lies above the line extends up to 
the upper bound of T^upper provided by flM|) ). Then, via simple algebraic manipulation 
with the r.h.s. of the estimate for T^upper, the assertion (b) of the lemma follows. 

Finally, since g'p{-) < b (as ffT^ ensures) then the graph of gp{d), d G [j,oo), lies 
below the line with slope b constructed by us. Hence, TZiower < T^{p), and the assertion 
(a) of the lemma follows via simple algebraic manipulation with the r.h.s. of the lower 
bound for T^iowcr provided by (j34l) . □ 
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Lemma 3. Let {pn, dn), n >0, be an orbit of (QP- Then, either {pn, dn) (0, 0) or there 
exists finite r > such that 

either < pr+i < Xdr or Xdr < Pr+i < 0. (35) 



Proof. Case 1: First, we consider the case when do and po have different signs. We shall 
conduct our arguments under the assumption that po > and do < 0; in the opposite 
case, i.e., when Po < and do > 0, the proof follows by exactly the same argument with 
the obvious change of inequality directions. 

The assumption po > and do < and the relation (fTSi) imply that: 

Pi < Po and pi > Xdo- (36) 

For this reason, if pi < then (1351) holds with r = and thus, the lemma is proved. 
Assume the contrary: pi > 0. Let us prove that lii < under this assumption. Since 

\d,-Vip,)\ ™^ \g,,ido)-g,mPi)\ ? Ido-Vip,)], (37) 

then \di — V{pi)\ < \do — T>{pi)\. The inequality (jSTj), the inequality 'E'(pi) < (that stems 
from flTTI) and the assumption pi > 0) and the inequality do < (the initial assumption) 
can hold simultaneously only if di < 0. This is the conclusion we aimed for. 

The argument of the above paragraph either finishes the proof of the lemma or implies 
that pi > and di < hold. In the latter case, the argument can be repeated with these 
inequalities in the place of the argument's assumption "po > and do < 0". Upon 
repeating this procedure, we either eventually find an r such that p^+i < or not. When 
the first alternative is the case, the last but one repetition of the argument ensured that 
Pr > 0. Combining this inequality with pr+i < and equation f|T5l) . we get that the 
second double inequality of fl35|) holds true for the found value of r, and thus the lemma's 
proof is finished. The second alternative is analyzed in the paragraph below. 

Let us assume that the argument of the last but one paragraph can be repeated 
infinitely. Since, after n-th repetition, we have that Pn > (otherwise, r would have been 
equal to n — 1) and dn < (this inequality is the final conclusion in each repetition) then, 
in virtue of (fT5l) . Pn+i < Pn and therefore, the sequence {pn}'^=o must have a non-negative 
limit; we denote it by p. The convergence Pn ^ P and the relation dn-i = \^^{pn —Pn-i) 
(that follows from (1151) ) ensure then that d^ — > 0. Therefore, (p, 0) is an equilibrium of 
([T]). Due to Theorem [H p = 0. This proves the lemma's assertion. 

Case 2: Now we consider the case where do and po have the same sign. We shall assume 
that Po > and do> 0; the opposite case can be treated by exactly the same argument. 
We observe initially that 

if for some n > 1, dn-i > and Pn-i > then p„ > Pn-i and d^ < dn-i- (38) 

The first conclusion in the implication (1381) stems directly from its assumption and from 
f|T5|) . For the second one, we argue as follows. If TZ{pn) 7^ ±oo then the following relations 
hold: 

JTBl I21t 12011 

Mnl = \gp„{dn-i) - gp,^ (n{pn))\ < |(i„_i - 7^(p„)| . (39) 
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From (jSHD, we obtain that \dn\ < \dn-i — 7^(p„)|. Next, from p„ > and (^Bi), we obtain 
that lZ{pn) > 0. The two inequahties just obtained and the assumption > of (1551) 
can hold altogether only if dn < dn-i- This completes the proof of (138|) in case 7l{pn) 7^ 
±oo. Let us analyse the opposite case. First, we note that the inequahty Pn > (which 
is the consequence of the assumption Pn-i > and our conclusion p„ > p„_i) and the 
property (l25ll altogether imply that TZ{pn) can only be +oo. This conclusion and property 
(IMl) ensure then that gp„{d) < for all d, and consequently, that dn = fi'p„(rfn-i) < 0. 
Since c?„_i > has been assumed then dn < rfn-i follows and fl38|) is established. 

Consider now an orbit {(Pn;C^n)}^o satisfying our assumption Po > and d^ > 0. 
Obviously, either rf„ > for all n, or we shall find m < oo for which (i^ < 0; there may 
be several such numbers, let m denote the smallest of them. For this m it holds that 
Pm > 0. This is ensured by recursive application of fl55]) for all n < m. Hence, the portion 
{{pn, dn)}'^=m of the Considered orbit can be treated by the argument of Case 1. By this 
argument, the lemma's assertion follows. It remains only to prove the lemma for the case 
when dn>^ for all n. We do this below. 

Assume (in addition to the assumption Pq>Q and (io > already made) that (i„ > 
for all n. The relation (1381) implies that in this case, is a non-decreasing sequence 
of non- negative numbers. If it is unbounded from above, then so is the sequence {TZ{pn)} 
in virtue of Lemma [2] (a). But since, by the very definition, \dn\ < 1 for all n, then there 
must exist k such that d^ < 7l{pk+i). For this k, 

rrm ( (ED. dk<Tl{pk+i) ||2^ 

dk+i = gp^^Mk) { < ^'Pfc+i (^(P'i+i)) = 0' ^(Pfc+i) is finite, 

[ < 0, since gp^+-^{-) < when TZ{pk+i) = +oo 

(we used above that lZ{pk+i) cannot be — oo when pk+i > 0) and therefore, dk+i < 0. 
This contradicts the current assumption that dn > for all n. Hence, {pn} must be 
bounded and thus, converges to a finite limit p. From this convergence we deduce that 
{pn, dn) — > (0, 0) in the same way as it has been done four paragraphs above. □ 

Lemma 4. k be the minimal positive integer for which < 1 -^"(i-^) ^^/^g condition 
b G (0, 1) guarantees the existence of such k). Let Ac be defined as the maximal positive 
real number satisfying the following inequalities: 

(0 Xck< Y'^^\''^\ (n)A.<aJ^, {iii)K<h. (40) 

Suppose A G (0, Ac]. Let {pn, dn),n > 0, be an orbit of Then, 
if ^ Pi < A(io, then either {pn, dn) (0, 0) or there exists a finite i > 1 such that 

(i) dj^i > Tl{pj), (li) 0< dj < bdj_i, (iii) pj+i > pj, for all j = !,■■■ , i, 
and (iv) di < TZ{pg+i)] 

if ^do < Pi < 0, then either {pn, dn) — > (0, 0) or there exists a finite i > 1 such that 

(i) dj_i < n{pj), (li) bdj_i < dj < 0, (lii) Pj+i < pj, for all j = !,■■■ , i, 
and (iv) dn > Tl{p£^i). 

^b, a and h employed here for the definition of k and Ac have been defined in ([T7| . ([30]) and ((32|) . 
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Proof. We shall conduct the proof under the assumption < pi < Xdo; in the case 
A(io < Pi < the argument is analogous and will not be presented. 

Step 1: We claim that the following inequality stems from lemma's assumptions: 

do > 7^(pl). (43) 

Indeed, by the very definition fl68l) . \do\ < 1 and hence the assumption < pi < Xdo 
implies that < pi < A. This, together with the assumption A < Ac and the constraint 
(140]- (iii)) yield that < pi < h and thus we can apply Lemma [2](b). It gives: 

7^(pl) < (44) 
aJ a 

On the other hand, the assumption pi < Xdo together with the constraint (HO] -(ii)) imply 
that do >-^^Pi- This inequality and (glD yield (03]). 
Step 2." We get: 

di = gpAdo) > (7^(pl)) = 0, (45) 

and 

\di\^= \gpMo)-9pAnPi))\ < b\do-n{pi)\<bdo, (46) 

where the last passage in (H6l) is valid because do > (this is one of the inequalities 
assumed in the beginning of the proof) and < 7?.(pi) < do (here, the first inequality is 
provided by (1231) and the assumption that pi > 0, while the second inequality is identical 
to (H5]) ): note also that this double inequality ensures that 7?.(pi) 7^ ±00 and hence, the 
last passage in both (H5l) and (H6|) is legitimate. 

Step 3: Since di > 0, as ensured by (H5l) . then (|T5|) yields that 



P2 > Pi. (47) 

This is the final conclusion of the third step. 

The relations (H5|) . (H6l) and (H7j) prove the vahdity of (i), (ii) and (iii) of (HT]) for j = 1. 

Now, if di < TZ{p2) then the lemma is proved with i = 2. Thus, we continue the proof 
assuming the contrary: di > 7l{p2). Obviously, this assumption implies that di > TZ{p2)- 
Taking this relation in the place of (H3!) and repeating the second and the third steps of 
the argument presented abov^ we deduce the relations < ^2 < bdi and ps > p2- These 
relations and the inequality di > Ti{j>2) imply that (i), (ii) and (iii) of (HTl) hold for j = 2. 

It is obvious that the argument of the above paragraph can be repeated for j > 2 
provided dj > 7l{pj+i. However, the repetition process cannot last forever, unless pi = 0. 



^Note that we do not need to repeat the argument's first step since its conclusion (P5| is now a direct 
consequence of our assumption. Hence, we will not employ Lemma[5][b) in the current and the consequent 
repetitions. Getting rid of the necessity for the use of this lemma is here, because we cannot guarantee 
that \pj\ < h, for j > 2, and therefore, we cannot ensure the validity of lemma's assumptions at j-th 
repetition for j > 2. 
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The reason for this is the following. After j consecutive repetitions, we would have 
that dj < WdQ (in virtue of fllTI-ii)) and pj+i > pi > (in virtue of the assumption 
Pi > and the property (HTI-iii)). Thus, the sequence {dj} decreases to zero, while the 
sequence {lZ{j)j)} possesses - in virtue of Lemma[2]^a) - the following property: TZ(j)j+i) > 
-^fPj+i > ^fPi- Consequently, if pi 7^ then there must be a finite i for which 
di < 7l{pe+i), and the lemma is proved. 

To finish the proof, we have to complete the argument of the above paragraph by 
analyzing the case pi = 0. Suppose in addition that do = 0. It then follows directly from 
([1]) that {pj,dj) = (0,0) Vj > 1. This conclusion completes the proof since one of the 
alternatives in the lemma's assertion is that (p„,(i„) — > (0,0). 

The last case to be considered is thus, pi = and do > 0. In this case, 5'pi(-) = fi'o(')- 
We note that (1171) and (El-d) ensure that 5-0(0) = ! - 2^{0)=0. Hence, we get: 



ifTgl „,_n do>0 and 1211 

di ^ 9pMo) goido) > goiO) = 0. (48) 

From f|T5l) and fHHl) . we get that p2 = Xdi > 0. This implies that < p2 < Xdi. With the 
latter double inequality in the place of < pi < Xdo we repeat the whole argument starting 
from the beginning of the proof and finishing at the end of the above paragraph. In the 
repetition, we shall not stumble upon the inconvenient possibility "p2 = 0" , since we have 
just shown that p2 > 0. The result is the following conclusion: either (p„, dn) — > (0, 0) or 
( 14T|) holds with j = 2, . . . , £ for some finite i > 2. As for the validity of (jH]-i, ii, iii) for 
j = 1, it has been already established (read the sentence after (H7I) ). This completes the 
proof of the lemma. □ 

Lemma 5. Suppose X G (0, AJ where Ac satisfies (J^- Let {pn,dn),n > 0, be an orbit of 
(Qy. Then, 

if < do < TZipi) then either {pn, dn) ~^ (0, 0) or there exists m > 1 such that 

(i) V{pi) < dj < 0, for j = 1, ■ ■ ■ , m, (ii) < Pj < Pj-i, for j = 2, ■ ■ ■ , m, ,^g- 
and (iii) Xd^ < Pm+i < 0; 

and if > do > T^ipi) then either [pn, dn) — > (0, 0) or there exists m > 1 such that 

(i) V{pi) > dj > 0, for j = 1, ■ ■ ■ , m, (zi) > Pj > Pj-i, for j = 2, ■ ■ ■ , m. 



and (iii) Xdm > Pm+i > 0. 



(50) 



Proof. We shall conduct the proof under the assumption < rfo < TZ{pi); in the case 
> rfo ^ T^iPi) the argument is analogous and hence will not be presented. 
First, we conclude that di < reasoning as follows. If TZipi) 7^ +00, then 



IJgji |2TJ, and the assumption do<7^(pl) 

di = gpAdo) < gpiC^ipi)) = 0. (5i) 

If TZipi) = +00 then in virtue of flM|) . gp^{d) < for all rf, and hence di = gp^{do) < 0. 
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Second, from fl2^ and the assumption Tl{pi) > we get that pi > and therefore, 
in virtue of Lemma [T], T>{pi) < 0. The latter and the assumption do > imply that 
do > 1^{pi)- We use this relation in the calculations below: 

1X51 l2U,do>»(pi) 1231 

and we conclude thus, that di > T'(pi). 

At the third step, we turn our attention to p2, but consider separately the cases P2 ^ 
and p2 > 0. 

Let us assume that P2 < 0. Our aim is to prove that in this case fHU]) is true for m = \. 
Since it has been proved in the previous two steps that T>{pi) < di <Q and since fHU]-(ii)) 
is void for m = 1, then our aim is achieved as soon as we prove that \di < P2 < 0. In this 
relation, the second inequality is exactly our current assumption, and the first inequahty 
is valid because p2 = pi + Xdi (the relation ( fTSl) ) and because pi > (proved at the second 
step above). The lemma is therefore proved in case P2 < 0. 

Let us now assume that p2 > 0. Our objective is to show that (l49]-(i)) and (l49] -(n)) 
are satisfied with j = 2. 

We start by establishing fH9] -(ii)) for j = 2. The first inequality is a direct consequence 
of our current assumption. The second inequality stems from the relation P2 = Pi + Arfi 
(the eq. (HM ) and the inequality di < that has been proved in the first step. 

We proceed by proving fH51-(i)) for j = 2. Using di < (proved in the first step) and 
P2 > (the current assumption) we get: 



d2 = gp2{di) < gp^{0) < gp^ i^—j < 0, 

that proves the second inequality in (l49] -(i)) for j = 2. Thus, to achieve our objective, we 
have only to show that 

d2 > V{pi). (53) 

We split the argument in two parts. First we prove (!53l) under the assumption that 
di > V{p2), and then under the assumption that di < V{j)2). 

Assume that di > V{p2). This yields the desired relation (!53l) as follows: 

iTgi ini, di>vip2) |23l 126J and o<p2<pi 

d2 = gpM > gp,{V{p2))=V{p2) > V{p,) (54) 

(note that the inequality < p2 < Pi used in the last passage is provided by fH9]-(ii)) for 
j = 2 that we have proved above). 

Suppose now that di < 'D{p2). Then 

,m GU, and the assumption di<2?(p2) HSJ s 

d2 = gp2{di) < 9p2\^yP2)) = T^{P2) (55) 

that ensures that d2 < V{p2). Next, we note that the assumption P2 > and the property 
(!27|) imply that T>{p2) < 0. Because of this inequality, the following two relations are 
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valid in the present case: di < V{p2) < and d2 < T^{p2) < 0. These relations and the 
inequahty 

\d, - V{p,)\ ® \g,,{d,) - g,, {V{p,))\ 9 \d, - V{p,)\ (56) 

can hold altogether only if di < d2. The latter and the relation T>{pi) < di imply (153|) . 

We can now conclude that if j92 > then (i) and (ii) of (H9l) hold for j = 2. Recall 
that in the first two steps of the proof we have already concluded that (l49] -(i)) is valid 
for j = 1. Combining these conclusions, we close our argument's third step: when P2 > 
then (i) and (ii) of (H9l) hold for m = 2. 

The third step of our argument, i.e., the proof that started right after fl52l) and finished 
above, can be repeated for P3,P4, etc. until we find m such that Pm+i < 0. The argument 
will then ensure that (jl9]-(iii)) holds for the value of m founded, and the previous recursion 
steps will ensure that fHU] -(i)) and fHm -(ii)) hold for this value of m. This finishes the proof 
of the lemma in case m is founded. If, to the contrary, such an m does not exist then 
the recursion ensures - via (H9]-(ii)) - that {pn}'^=i is a monotone non-increasing sequence 
bounded below by 0. Let p denote its limit. Then, from the equation (fT5|) we get that 
dn = (A)^^ {pn —Pn~i), and therefore, (i„ — > as n oo. But, in virtue of Theorem [1], 
the — > oo limit of (p„, dn) can only be (0, 0). This completes the proof of the lemma. □ 

Proof of Thm. |4l. In virtue of Lemma[3l in order to establish the theorem, it is sufficient 
to prove that {pn,dn) (0,0), in case when either \do > pi > or \do < pi < 0. We 
shall assume that 

Xdo >pi>0 (57) 

for the rest of the proof. In the case Arfg < Pi < 0, the proof is the same up to obvious 
changes of the inequality directions. 

The assumption flSTl) allows us to apply Lemma HI Let i be the integer provided by 
it (if i does not exist then (p„, dn) (0, 0) in accordance to the lemma and thus, the 
theorem is proved). 

First, let us prove that 

1 Ja{l -b) 

Pi < 2 1 ^0- (58) 

For the proof, we shall need the fact that pe > 0. It follows from the assumption pi > 
via a recurrent application of the relation (jUJ-iii) (that is valid due to Lemma H]). The 
continuation of our argument depends on whether (i — l) > k or (i—l) < k, where, recall, 
k has been defined in Lemma H] as the minimal positive integer for which 6^ < 1 . 
note that k depends solely on the parameters of the studied dynamical system, and hence 
can be applied to any its orbit. 

If (£ — 1) > A; we argue as follows. The inequality pe > just proved and the property 
yield the inequality pe < JTZ{pi), from which we get flSSl) via the following chain of 
estimates: 

Pi < J7l{pe) < Jde^i (due to (HlJ-i) of Lemma H for j = £) 

< Jb^~^do (from (l4T] -ii) of Lemma H] applied for j = l, 1) 

< Jb'^do (because (^ — 1) > A; and b < 1) 
— I -^"(1-'') (^lyy ^i^Q definition of k). 
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If {£—!)< k, then fl55]) follows from the following chain of inequalities: 





1 \ 7 1 

Pi + Adi + ■ ■ 


1 \ 7 

■ + Xdi_i 


(by recurrent apphcation of lllop j 


< 


Xdo + Xdi + ■ 


■ ■ + Xdi-i 


(using the first inequality of (1571)) 


< 


iXdo 




(using (HlJ-ii) of Lemma H] for j = 1, ■ ■ ■ ,i — 








and taking into account that 6 < 1) 


< 


kXdo 




(because (£ — 1) < /c in the considered case) 


< 


iJa(l-fe) , 
2 b "0 




(by the condition fHO]-(i)) on A). 



1, 



Second, from the inequality pi> just established and Lemma [T] we get that 

Vip,) > -—^p,. (59) 
a{l — 0) 

Third, we combine the conclusions fl58|) and fl59|) to deduce that 

V{pi) > -^do. (60) 

The fourth step is based on the double inequality < dg < 7l{pe^i); its first part is 
ensured by (HBii) with j = i and the second part by (jUJ-iv). We take this inequality 
in the place of the assumption of Lemma [5] and derive then from the lemma that either 
(pn, dn) — >■ (0, 0) and therefore the theorem is proved, or 

there exists m > 1 such that V{p£^i) < dj < 0, for j = i + 1, ■ ■ ■ , i + m, (61) 

and Xde+m < Pi+m+i < 0. (62) 

Combining (l60l) and (|6T|) . we get that 



-do < di+j < 0, for j = 1, . . . , m. (63) 



The inequalities ( l63l) just derived, the inequalities (HTT-ii) for j = and the 

inequality b < 1 lead altogether to the following conclusion: 

in the block di, . . . , di, di+i, . . . , di+m of the sequence {(i„}^0' 

the values of the first i{£ > 1) members are between and do (64) 

while the values of the last m(m > 1) members are between — ^do and 0. 

The argument that started at ( 1571) and finished at ( l64l) can be repeated with the 
inequahty (l62l) in the place of (!57|) . The respective conclusion is that the block di, . . . , c/^+m 
is followed by another block - of the size 2 at least - that consists of two non-empty parts 
such that: the members of the first part are all between dg^m and 0, and the members of 
the second part are all between and —^dg^m- This implies in virtue of the inequality 
M^+ml ^ \do (that stems from (Elj)), that the absolute value of each members of the block 
does not exceed | do- 
lt is obvious that the argument can be repeated yielding at the i-th step the conclusion 
that the members of the corresponding block of the sequence di,d2,... do not exceed 
(I)* do, in modulus. From this, dn follows. 
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To complete the proof, it is only left to show that p„ ^ 0. Suppose that this is not 
the case. Then, there exists e > and a sequence of integers {^^j}jgN such that rij t oo 
and \pnA > £ for each j. We suppose, without loss of generality that p„ . > e. Thus: 



dn, = gp„^{dn,_J < Qeidn,-^)- (65) 

The inequality ( l65l) . the continuity of 5'p(-)0 and our conclusion dn imply altogether 
that < ge{0). But, on the other hand, ( 12T1) and (l22l) ensure that > ge{0). Hence, by 
contradiction, we conclude that p„ 0. 

We thus have proved that dn ^ and p„ — for arbitrary po and c/q. This is the 
theorem's assertion. □ 



3 Application 

3.1 Our HI A model for a single risky asset market 

The bridge that leads from our results to their applications is grounded on the Hetero- 
geneous Interacting Agent Model to be presented now. We shall use the abbreviation 
"HIAM" . When we need to distinguish it from other HIAMs, we shall refer to it as "our 
HIAM" . 

The model's time set is {0, 1, 2, . . .}. The model's parameter denoted by K and called 
the number of model's agents is an arbitrary natural number. There are abstract agents 
in the model; they are labeled by the numbers 1,2, . . . , K. The other model's parameter, 
denoted by a and called the proportion of speculators, is an arbitrary real number in the 
interval (0, 1). It plays the following role in the model: the agents numbered 1,2,..., aK 
are called speculators, and the rest of the agent population, i.e., the the agents numbered 
aK + 1, aK + 2, . . . , K, are called fundamentalists. The choice of these names will be 
justified a few paragraphs below. 

To each time n > 1 and to each agent labeled k, the model associates a random 
variable denoted by dn{k). Its definition is as follows (the quantities involved in (!66|) and 
(|67ll will be specified below; in particular, note that Vn{k) will be a random variable - this 
is the source of the randomness of dn{k)): 

d(k) = -l ' ^' " lPn-iv + Vn{k))\>0 if k = 1,2,..., aK, 

' 1 ;f 7^ Sp^ _ _|_ ^^(^/j^^j < Q i_g_^ if is a speculator. 



+1, 


if 


Jdn-1 — 


-1, 


if 


Jdn-1 — 


+1, 


if 


- hn - 


-1, 


if 


- [Pn - 



d ()^) ^ ) ' yn yu+Vn{k))\>0 li k = aK + 1, aK + 2, . . . , K , 

\ —1, if — [pn — {v + Vn{k))\ < i.e., if A; is a fundamentalist. 

dn{k) is called the decision of agent k at time n; when dn{k) is +1 we say that "the agent 
wishes to buy an asset share", and when it is —1 we say that "he wishes to sell an asset 
share" This "buy/sell" interpretation helps to understand the link of the model with 



is continuous because it is difFerentiable - see (|18l) . 
^^In our HIAM, there is no an equivalent of the market clearing condition. Therefore, the one who 
wishes to buy/sell will not always accomplish this wish. By this reason, we why we say "wishes to 
buy /sell" and not simply "buys/sells". 
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the asset market and the names that we shall attribute to the model's parameters and 
variables introduced below. 

Let us specify the quantities involved in (166|) and (167|) . 

First, 

1 ^ 

dn:=j^J2dn{k), n = l,2,..., (68) 

k=l 

and in accordance with this definition and with the interpretation of each dn(k), the real 
world analogue of dn is the population relative excess demand for the asset at time n; the 
name for the analogue will be the name for dn in the model. However, when possible, we 
shall shorten this name to just the excess demand. Observe that (1^ defines dn for n > 1. 
As for do, it is one of the model's parameters called the initial excess demand; its value 
range is [—1, 1]. The value of Jo as well as of po to be defined below, must be specified in 
order to "switch on" the model's evolution. 

Second, Pn is the model's quantity that corresponds to the price of one share of the 
asset in the market mimicked by the modelQ We call Pn the asset price at time n, or 
simply price. The value range of pn is M for each time n = 0, 1, 2, . . .. 

The model's price evolves in time. The evolution rule is as follows: po is one of the 
model's parameters (it is called the price initial value), and all others p„'s are determined 
by the recurrent relation 

Pn+i =Pn + A4, n = 0, 1, 2, . . . , (69) 

where A is a positive real model's parameter called the feedback of the excess demand on 
price increment. The recurrent relation (1691) called price update rule mirrors the demand- 
and-supply law of real markets. Indeed, when dn is positive [resp., negative] - which 
means that there are more buyers than sellers [resp., sellers than buyers] - the rule rises 
[resp., lowers] the price. 

Third, we specify and interpret v and f„(fc)'s. w is a real value which is one of the 
model's parameters, and {f„(fc), n = 1,2,..., k = 1, 2, ... , K} are independent among 
themselves and of all other random variables present in the model, with a common distri- 
bution function that will be denoted by $; this function is another model's parameter. We 
assume that $ satisfies the the conditions ([3]). The value v + Vn{k) is called the evaluation 
of the asset's fundamental value by agent k at time n, or, simply, individual evaluation. 
The name refiects the quantity that we want to mimic by v + Vn{k) in our model, namely: 
what the agent k thinks the asset worths at time n, or as in traditional economical terms, 
the asset's fundamental value evaluated by agent k at time n. 

In the real world, asset's fundamental value is usually calculated from the official 
balance reports of the company that issued the asset. However, an individual evaluation 
of this value may not coincide with the officially calculated one because an individual may 
have privileged information (about the future company development, say), or may use his 
own way to deduce the fundamental value from the report data. This situation is mirrored 

""^^ Actually, it corresponds to the to the logarithm of the price of one share of asset. By using the 
logarithm, we allow the model's counterpart of price be negative. This simplifies mathematical treatment 
of the model. Nevertheless, the name for p„ in the model will be "price" rather than "logarithm of price" . 
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in our model in the following way: v corresponds to the official asset's fundamental value 
that may be simply the average taken over all individual evaluations, and differences 
between this value and individual evaluations are modeled by the random variables f 's 
with the distribution function $. Accordingly, the name for v in our model is the market 
fundamental value of the asset, and the name for $ is the distribution of deviation around 

V of individual evaluations of the fundamental value of the asset. 

Up to now, we have introduced and interpreted all the quantities involved in the 
decision rule (|67|) . This allows us to make the following observation: an agent that "acts" 
due to this rule, will wish to buy one asset share if he thinks that the current price (i.e., 
Pn) is less than what the asset share worths in accordance to his current evaluation (i.e., 

V + Vn{k))] otherwise, he will wish to sell one asset share. The one who behaves in this 
way in the real market, is called fundamentalist. This explains why we use this name for 
our model agents that obey the rule (1671) . 

Finally, we explain the role of J. This model's parameter is a positive real number 
that mirrors the social susceptibility, or simply susceptibility of market agents. It models 
"susceptibility" because its use in the rule (1^ makes it to correspond to the strength 
with which the population excess demand influences an individual to align his decision 
in accordance with this excess. Now observe that the population excess demand is the 
average of decisions taken over the whole population and hence the alignment means the 
coincidence of an individual decision with that of the majority of the population. Hence, 
the "susceptibility" mirrored by J has a "social" aspect of an individual behavior. 

The parameter J is also interpreted as the traders' speculative trend. The justification 
for this is as follows. In our model, the population excess demand is proportional to the 
price trend (by the price trend we mean pn+i — Pn, and the proportionality is yielded by 
(!69|) ). Consequently, the agent who ahgns his decision with the population excess demand 
acts as an asset market speculator (i.e., an individual that wishes to buy [resp., sell] when 
the trend indicates the price will rise [resp., fall]). This fact explains why the model's 
agents that "use" the decision rule (IMi) are called speculators. However, these agents 
are not pure speculators because the rule fl66l) mixes the speculative behavior and the 
fundamentalist behavior (explained two paragraphs above). In this mixture, the weight 
of the speculative component is given by J. This role of J suggests its interpretation as 
the market traders' speculative trend. 

3.2 Why and how the studied dynamical system mimics evolu- 
tion of asset price and excess demand in asset markets 

Corollary [1] will answer the question posed in this section title. It grounds upon Theorem[5] 
and Lemma [6l both to be stated below. We could state that the convergence affirmed 
in the theorem holds true because the Strong Law of Large Numbers ensures that dn 
converges to a constant, as the number of agents increases. However, it is impossible to 
construct a rigorous proof based solely on this law. Some adequate estimates for deviations 
are necessary. The idea behind them is simple but the execution requires heavy notations 
and cumbersome calculations. An interested reader may find all this in p.Oj. As for the 
proof of the lemma, it can be conducted easily by induction in time n, and hence will not 
be presented here. 
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Theorem 5 (On convergence of trajectories of our HIAM to orbits of a ([TJ-like 
dynamical system). 

Choose a distribution function $ satisfying the conditions ([^^ and real numbers a G (0, 1), 
X> 0, V eR, J > 0, po and do e [-1, 1]. 

For each K = 1,2,..., and the entities chosen above construct the HIA model as 
defined in Section \3.1\ and denote the price and the excess demand at time n of the con- 
structed HIAM by p^'^ and d^'^ , respectively^ 

Next, forn > 1, define the sequence of pairs of real numbers {{pn , d^) , n = 0, 1, 2 . . .} 
so that do = do and for n > 1, 

( Pn= Pn~l + Xdn-1 

\ c?, = «[l-2$(p„-t;-Jc/„_i)] +(!-«) [1- 2$ ^ > 
Then, for any t eN and any e > 0, it holds that 

lim Pfldf ^ - rf„| < e, and\p'^'^ -pn\ <e, for all n = 0,1,2, ... ,t) = 1. 



It can be easily noted that (!70|) . the dynamical system that figures in Thm. |5l is slightly 
different from (Q, the dynamical system that we study in Section [2l Nevertheless, as the 
following lemma reveals, trajectories of these systems have a very simple and rigid link. 
This link allows one to "transport" any result about ([T]) into the framework of (170 p . 

Lemma 6 (on the link between the dynamical systems ([T]) and ( 1701) ). 

Suppose the dynamical systems fW\ l and (J\) have the same values for the parameters a, 
X, J, and the same distribution function $ and let the value of the parameter v of ( fTOl ) be 
arbitrary. Let us add the superscript * to Pn and dn in the system (70) so that they might 



be distinguished from pn and dn in the system (Qp. Suppose that p^ and dg, the initial 
conditions for an orbit of (7(^ , relate to po and do, the initial conditions for an orbit of 
(Qlj, as follows: Po = Po — v, and do = dQ. 

Then, for every n>l, it holds that Pn = Pn — v and dn = • 

Thm. [5] and Lemma [6] yield straightforwardly the following result. 

Corollary 1 (why and how the dynamical system ([T]) mimics evolution of asset price 
and excess demand). 

Let Po E M., do E [—1, 1], a G (0, 1), J > 0, A > 0, w G M, and a probability distribution 
function $ satisfying the conditions ^ correspond to characteristics of a single risky asset 
market as specified in the Introduction. 

Consider the dynamical system (Q]) determined by a, J, A and $. For arbitrarely fixed 
time T, let {{pn, dn) n = 0,1, . . . , T} denote the orbit of the system starting from {po, do). 

Then, the trajectory on the time interval [0, T] of the difference between the price of the 
market asset and the asset's fundamental value lies closely to {pn, n = 0,1, . . . , T}, while 



^•^The purpose of these modifications is to emphasize that the price and the demand in the HIA model 
constructed in Section [3?T1 are stochastic quantities and that they depend on K, the model's population 
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the trajectory on the time interval [0, T] of the excess demand for the asset lies closely 
to {dn, n = 0,1, . . . ,T}. Moreover, the larger the market agent population size the closer 
are the trajectories to their respective approximations. 

3.3 Revealing and explaining market price and population ex- 
cess demand dynamics properties 

In the present section, we employ Corollary [T] to interpret the results of Section [2] in terms 
of a real world asset market. For the interpretations to be valid, it is necessary that T 
and the market agent population size (both from Corollary [1]) be large. We take it for 
granted that both are as large as necessary, when necessary. 

Before we proceed, it must be noted that assumptions of the results of Section [2] do 
not concern the whole distribution function $ but rather solely $'(0) (that is, the value 
of the derivative of $ at 0). This is an interesting virtue of our results. However, when 
one wants to apply them to explain real world market properties this virtue becomes an 
obstacle since there is no a generic interpretation for $'(0) in terms of real world markets 
(that is, an interpretation that would suit any $). To overcome this obstacle, we shall 
accept the following assumption: 

$ is a Normal distribution function with zero mean. (71) 

The point here is that if o"^ denotes the variance of $ then this assumption implies that 
0"^ = (-\/27r<I>'(0)) which, in turn, allows us for the following interpretation: 

($'(0))^^ corresponds to the heterogeneity of the distribution of individual 
evaluations of the asset's fundamental value, in the sense that larger value of (72) 
($'(0))~^ corresponds to higher dispersion of evaluations over the population. 

The interpretation fl72l) is extremely convenient in re-phrasing the results of Section [2] 
in terms of real world asset markets. This is a strong motivation for accepting the as- 
sumption ( ITTI) . There are other motivations. One of them stems from our belief that 
an individual is influenced by diverse factors and information streams, when he makes 
up his mind in respect to the fundamental value of an asset. Consequently, due to the 
Central Limit Theorem, individual evaluations of this value should be distributed over 
a population in accordance to a Normal Law. Another motivation comes from the fact 
that the Normal Law would be the most convenient and robust theoretical model for a 
populational distribution obtained by sampling from a real world population. In other 
words, the Normal Law would likely be chosen for $ when fitting our HIAM to a real 
world socio-economic process. 

Corollary 2 (of Thm. [1] obtained with the help of Cor. [1]). 

(a) The market state in which 

the market asset price coincides with the asset's fundamental value . 
and the excess demand for the asset is zero 
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is the only possible equilibrium state for the asset price and the excess demand of an asset 
market. 

This means, in particular, that if the asset price and the excess demand converge to 
some value as time goes on, then their limit values are respectively, the asset's fundamental 
value and 0. 

(b) (valid under additional assumption (|7T]) ) Whether the state (7^ is a locally stable 



or an unstable equilibrium is determined by a relation between two expressions that we 
present below and denote by u and w: 

u := 2 X (proportion of speculators) x (traders' speculative trendjx 
X ( the heterogeneity of the individual evaluations 

of the asset's fundamental value)~^ 
and (74) 
w := 2 X (the feedback of the excess demand on price incrementjx 
X ( the heterogeneity of the individual evaluations 

of the asset's fundamental value)~^. 

The determinant relation is as follows (it is illustrated in FigureU\): if 

w < 2u + 2 and u < 1 

(together with u > 0,w > which are implicit from [74\ ) and the natural 
constraints that impose that each entity in the definition [74\ ) is positive) 

then ([73y is locally stable, otherwise, it is unstable. 



Corollary 3 (of Thm. H] obtained with the help of Cor. [T])). 

There exists a strictly positive threshold Ag such that when the feedback of the excess 
demand on the price increment is smaller than Xq then the local equilibrium described in 
Corollary d turns to be global, that means that the long time limit values for the asset 
price and the excess demand for this asset are, respectively, the asset's fundamental value 
and 0, whatever the initial values are. 

Corollary 4 (of Thm. [3] and Conjecture [2] obtained with the help of Cor. [T] and valid 
under additional assumption flTTl) . 

Suppose that, as a result of a smooth change of the parameters' values, an asset market 
is transferred from the parameter values region in which its equilibrium is stable to the 
region in which its equilibrium is unstable. Suppose that the change is such that 

(^ij (a) the traders' speculative trend (i.e., the value of parameter J) increases, or 

(b) the proportion of speculators (i.e., the value of parameter a) increases, or 

(c) the heterogeneity of the distribution of the individual evaluations of the asset's 
fundamental value diminishes (i.e., the value of^'{0) increases), or 

(d) any combination of (a)-(c). 

(a) the feedback of the excess demand on the price increment (i.e., the parameter X) is 
kept fixed, and its value is sufficiently small so that X/{aJ) < 2 for all values of a 
and J during the process of change. 
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Then, in the transfer course, there will be a (maybe short) interval of time during which 
the market asset price and excess demand exhibit regular oscillation with non damped 
amplitude. 

We would like to close the presentation with several comments: 

(a) The last statement of Corollary [2] is important because a real world asset market 
always experiences endogenous and/or exogenous shocks, and therefore the stability prop- 
erties of the equilibrium (1731) become an important issue. 

(b) In real world markets the price update increments are small, hence Corollary [3] sug- 
gests that, in the real world, when the state fl73l) is stable, it is actually globally stable. 

(c) One can conclude from Theorem [3] and Conjecture [2] that for a specific combination 
of market parameters' values, there will appear oscillations of asset price and excess de- 
mand for that asset. This conclusion has limited application since the precise estimation 
of market's parameters is a very difficult task. Contrasting, the conclusion presented in 
Corollary |4] has a practical value because it describes a qualitative aspect of market be- 
havior. From this description, one indeed can derive useful results, like the following one: 
Let one know that the corollary's assumption (i), (ii) hold true and suppose one observes 
that a market exhibits oscillations of the asset price and the excess demand. Then, even 
though the oscillations might have disappeared, one can affirm that the market is in an 
unstable state. 
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